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Abstract. Using the HEGRA system of imaging atmo- 
spheric Cherenkov telescopes, a region of the Galactic 
plane (-10° < b < 5°, 38° < I < 43°) was surveyed for 
TeV gamma-ray emission, both from point sources and of 
diffuse nature. The region covered includes 15 known pul- 
sars, 6 known supernova remnants (SNR) and one uniden- 
tified EGRET source. No evidence for emission from point 
sources was detected; upper limits are typically below 
0.1 Crab units for the flux above 1 TeV. For the diffuse 
gamma-ray flux from the Galactic plane, an upper limit 
of 6.1 • 10 -15 ph cm~ 2 s _1 sr _1 MeV -1 was derived under 
the assumption that the spatial distribution measured by 
the EGRET instrument extends to the TeV regime. This 
upper flux limit is a factor of about 1.5 larger than the 
flux expected from the ensemble of gamma-ray unresolved 
Galactic cosmic ray sources. 
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1. Introduction 

Systems of imaging atmospheric Cherenkov telescopes, 
such as the HEGRA stereoscopic telescope system (Daum 
et al. 1997, Konopelko et al. 1999a), allow to reconstruct 
the directions of air showers over the full field of view - 
with a radius of about 2° - and can therefore be used 
for sky surveys (Puhlhofer et al. 1999). Here, we report 
on a survey of a rectangular patch of the sky of roughly 
80 deg 2 , centered on the Galactic plane at longitude 40°. 
The motivation for this survey was twofold: 

— Search for diffuse gamma-ray emission from the Galac- 
tic plane 

— Search for gamma-ray point sources 

Diffuse emission from the Galactic plane results from the 
interactions of charged cosmic rays with interstellar gas 
confined to the plane or with photons. Diffuse emission in 
the energy range from tens of MeV to tens of GeV has 
been studied intensively by SAS 2 (Fichtel et al. 1975; 
Hartmann et al. 1979), COS B (Mayer-Hasselwander et 
al. 1980, 1982) and EGRET (Hunter et al. 1997). The 
basic features can be modeled assuming ir° decay as the 
dominant mechanism, with the gamma-ray emission pro- 
portional to the product of the gas column density and 
the cosmic-ray density (see, e.g., Hunter et al. 1997, Fich- 
tel et al. 1975, Strong et al. 1988, Bloemen 1989, Bertsch 
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et al. 1993). The distribution of cosmic rays is assumed to 
follow the matter density with a characteristic correlation 
scale of 1.5 to 2 kpc (Hunter et al. 1997). In addition to 
diffuse Galactic gamma rays, there is also a small extra- 
galactic component, which should be fairly isotropic (Fich- 
tel et al. 1978; Sreekumar et al. 1998). Above 1 GeV, data 
show an excess in gamma-ray flux over model predictions 
(Hunter et al. 1997; see however Aharonian and Atoyan 
2000). At these energies, contributions from inverse Comp- 
ton scattering of electrons start to become relevant (see, 
e.g., Porter & Protheroe 1997). In response, revised mod- 
els speculate that the local measurements of the elec- 
tron flux may not be representative for the entire Galaxy 
(Porter & Protheroe 1997; Pohl & Esposito 1998). Elec- 
tron propagation is limited by radiative losses, and the lo- 
cal electron spectra arc strongly influenced by the history 
of sources in the local neighborhood (Aharonian, Atoyan 
& Volk 1995; Pohl & Esposito 1998). In case that the 
solar system is in an "electron void" , diffuse gamma-ray 
emission at high energies could be an order of magnitude 
above predictions based on local electron spectra (Pohl & 
Esposito 1998). Detailed models of the full spectrum of 
diffuse gamma-ray emission from the Galaxy (Moskalenko 
& Strong 2000; Strong, Moskalenko & Reimer 1999, 2000) 
also favour a harder electron spectrum. In the Galactic 
plane another 'diffuse' gamma-ray flux component arises 
from the hard energy spectrum of those Galactic CRs 
that arc still confined in the ensemble of their unresolved 
sources. Assuming these to be SNRs, Berezhko & Volk 
(2000) showed that their spatially averaged contribution 
to the diffuse gamma-ray flux at 1 TeV should exceed the 
model predictions of Hunter et al. (1997) by almost an or- 
der of magnitude. It is therefore of significant interest to 
search for the extension of the diffuse emission from the 
Galactic plane at higher energies. Upper limits on diffuse 
gamma-ray emission at TeV energies have been reported 
by Reynolds et al. (1993), and LeBohec et al. (2000), at a 
level of a few 10 -3 of the cosmic-ray flux. The most strin- 
gent limits at higher energies, above 100 TeV, were given 
by Borione et al. (1997), and constrain the diffuse flux to 
less than 3 ■ 10~ 4 of the cosmic-ray flux. 

The Galactic plane is also a region rich in potential 
gamma-ray point sources. For obvious reasons, supernova 
remnants as well as pulsar driven nebulae cluster along 
the Galactic plane. Both types of objects are almost cer- 
tainly accelerators of cosmic rays and emitters of high- 
energy gamma radiation. Theoretical models predict that 
typical gamma-ray fluxes from the majority of these ob- 
jects are below the detection threshold of the current gen- 
eration of instruments (see, e.g., Aharonian et al. 1997; 
Drury et al. 1994). However, both the lack of knowledge 
of the individual source parameters as well as approxi- 
mations used in the modeling result in large uncertain- 
ties in the predictions for individual objects, by an order 
of magnitude or more. In addition to pulsars and super- 
nova remnants, many unidentified EGRET sources lie in 



the Galactic plane. Given the density of source objects, a 
survey provides an efficient way to search for gamma-ray 
emission. 

The selection of the region to be surveyed was governed 
both by astrophysical and by practical considerations. Us- 
ing the EGRET observations of the diffuse gamma-ray flux 
(Hunter et al. 1997) as a guideline, a region within the 
inner Galaxy < 40°) and close to the Galactic equa- 
tor should show the strongest emission, corresponding to 
the gas column density. Similarly, the density of super- 
nova remnants (Green 1998) and pulsars (Tayler 1993) is 
highest in this region. However, from the location of the 
HEGRA telescope system, at 28°45' N, the Galactic cen- 
tre region can only be observed at rather unfavourable 
zenith angles of 60° or more, and the observation time 
is rather limited. At such large zenith angles, the energy 
threshold of the telescopes is increased from 500 GeV for 
vertical showers to 5 TeV at 60° (Konopelko et al. 1999b). 
Best sensitivity is obtained for observations at zenith an- 
gles below 30°, favouring regions at larger Galactic longi- 
tude. Another potential problem for Cherenkov observa- 
tions is the large and spatially varying background light 
from regions of the Galactic plane, which may result in 
non-uniform sensitivity across the field of view of the 
Cherenkov telescopes. 

As a result of these considerations, a rather dark re- 
gion at about 40° longitude was chosen for the survey, 
with an extension in Galactic longitude of 5° and in lati- 
tude of 15° (Fig. ||). The relatively large range in latitude 
was chosen in order to cover the full range expected for 
diffuse gamma-ray emission, taking into account that the 
inverse Compton mechanisms may result in a wider lati- 
tude distribution than observed in the EGRET data, and 
to include additional background regions at large latitude. 
The observation region includes part of the "Sagittarius 
Arm" , one of the spiral arms of the Galaxy. In addition, 
it hosts 15 known pulsars (Table ^), 6 supernova rem- 
nants (Table ||) , and the unidentified EGRET source 3EG 
J1903+0550 (Hartman et al. 1999). There are no bright 
stars in the survey region, which might seriously impact 
the observations; the brightest star is of m v = 4.3™ aff , and 
there are 10 stars brighter than 6 mag . (A 4.3 mQff star in- 
creases the DC current in a camera pixel from 0.8 to 10 
//A, resulting in an increased noise of 1.5 ph.e. rms). 

2. The HEGRA I ACT system 

The HEGRA stereoscopic system (Daum et al. 1997, 
Konopelko et al. 1999a) of imaging atmospheric 
Cherenkov telescopes (IACTs) is located on the Canary 
Island of La Palma, on the site of the Observatorio del 
Roque de los Muchachos, at 28°45' N, 17°53' W, 2200 m 
a.s.l. The stereoscopic telescope system comprises of five 
telescopes (CT2-CT6). One additional telescope (CT1) is 
operated in stand-alone mode. The system telescopes are 
arranged on the corners and in the centre of a square with 
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Fig. 1. Mean scaled width w of Cherenkov images in 
events where at least four telescopes triggered, for (•) 
gamma-rays from the Crab reference sample, and for (o) 
cosmic-ray showers. 



100 m side length. Each of them is equipped with tes- 
salated 8.5 m 2 mirrors of 5 m focal length, and a camera 
with 271 photomultiplier pixels in the focal plane. The 
field of view of each camera has 4.3° diameter; the pixel 
diameter corresponds to 0.25°. The telescope system is 
triggered when in at least two cameras two neighboring 
pixels show a signal above 6-8 photoelectrons. Signals from 
the cameras are recorded using a 120 MHz Flash-ADC 
system, which is read out after a trigger. Details of the 
camera hardware and of the trigger system are given in 
Hermann (1995) and Bulian et al. (1998). The pointing of 
the telescopes is known better than 1 arcmin. (Piihlhofer 
et al. 1997). On the basis of the stereoscopic analysis of 
Cherenkov images, shower directions can be reconstructed 
with a precision of 0.1°, and shower energies with a reso- 
lution of 20% or better (Daum et al. 1997; Aharonian et 
al. 1999a, b). For vertical incidence of showers, the energy 
threshold is 500 GeV for gamma-rays, and increases to 0.9 
TeV at 30° and 1.8 TeV at 45 (Konopelko et al. 1999b). 
Cosmic ray showers are suppressed exploiting the width 
of Cherenkov images. A "mean scaled width" w is defined 
by scaling the observed widths to the expected width for 
gamma-ray images, which depends on the intensity of the 
images and the distance to the shower core, and averaging 
over telescopes. Gamma-rays cause a peak in w at 1, with 
a Gaussian width of about 0.1 to 0.12 (Fig. |l|). Nucleonic 
showers have larger w values, peaking around 1.7. While 



more sophisticated identification schemes (e.g., Daum et 
al. 1997) can reach slightly better sensitivity, the default 
(and most stable) analysis schemes are based on cuts in 
w (w < 1.1. ..1.3), combined with an angular cut relative 
to the source of about 0.1° 



3. The data set 

The data used in this survey were taken during 37 days in 
June, July and August 1999 with the complete 5-telescope 
system. 

As illustrated in Fig. |^, the total observation time of 
88 h was distributed between three scans in Galactic lat- 
itude, centered at Galactic longitude of 39°, 40.5° and 
42°, chosen to guarantee an overlap of the effective fields 
of view for the three scans (the scan positions given here 
and below refer to the center of the field of view). Each 
scan was conducted in 1° steps in latitude between —4° 
and +4°, with an additional control region at —9°. The 
typical time per scan point was about 15 min per day, 
and the scans were repeated on several days. The points 
at latitude 0° received twice the exposure, once covered 
going from +4° to 0° and once going from —4° to 0°. A 
quality selection of the data sets was based primarily on 
the average trigger rate of the telescope system. A fraction 
of the data set suffered from extinction of Cherenkov light 
due to Sahara dust in the atmosphere (the so-called Cal- 
ima, a well-known phenomenon at the Canary Islands). 
The remaining data set is quite uniform in detection rate 
and not affected by Calima. In total it encompasses 41.7 
h, 18.6 h for the latitude scans at 39° longitude, 14.6 h at 
40.5° longitude and 8.5 h at 42° longitude. The data cover 
zenith angles between 21° and 35°; the zenith angle of a 
run and the Galactic latitude are somewhat correlated, 
with data at 4° latitude covering zenith angles between 
20° and 28°, compared to 28° to 35 for the runs at -9° 
latitude. In total, 1.4 • 10 6 Cherenkov events are used. 

The analysis of Cherenkov images could potentially 
suffer from variations in the sky brightness over the scan 
region. Since the readout electronics of the telescopes is 
AC coupled, a star illuminating a pixel will not cause 
baseline shifts, but it will still result in increased noise 
in that pixel. As a measure of the influence of night sky 
background light, the baseline noise of the photomultiplier 
signal was determined. The RMS of the baseline noise of 
the individual PMTs is quite homogeneous over the re- 
gion and amounts in average to about 1 ADC count. The 
brightest stars in the region increase the RMS to 1.5 ADC 
counts; this is uncritical for further image analysis. For ref- 
erence and comparison, observations of the Crab Nebula 
in the winters 1997/1998, 1998/1999 and 1999/2000 were 
used, and were subject to identical selection criteria. In 
total, 41.7 h of observations at zenith angles between 18° 
and 32° were used, with a roughly uniform coverage of the 
zenith angle range. 
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Fig. 2. Borders of the survey region (solid lines), and ob- 
jects in this region. The borders are shown for an effective 
field of view of 1.7° radius, corresponding to the radius 
where the detection efficiency is half of that compared to 
the value on the optical axis of the telescopes. Individual 
scan points are indicated as dashed circles. Pulsars are 
shown as black circles, supernova remnants as shaded re- 
gions, and stars brighter than 6 maff as crosses. The num- 
bers of pulsars and SNR correspond to Tables ^ and ||. 
The light gray area indicates the position of the EGRET 
source 3EG J1903+0550 and the la error box. 



The last telescope integrated into the system, CT2, 
one of the corner telescopes, was only used for part of the 
Crab reference data sets. Therefore, all analyses discussed 
in the following were performed both using only the four 
telescopes, as well as the full set of five telescopes. Because 
of the good agreement between the simulations and the 



Table 1. Known pulsars in the survey region. Columns 
include the Galactic longitude I and latitude b, the period 
P and the period derivative P. The dispersion measure 
(DM) can be used to estimate the distance to the pulsar, 
assuming a uniform density of thermal electrons in the 
Galaxy, in the order of 0.1-0.01 cm -3 and dividing the 
DM by this number. Data taken from Tayler (1993). 





Name PSR 
J(2000) 


in 


b[°] 


p N 


P 
[lO" 15 ] 


DM 
[cm _3 pc] 


1 


1909+0254 


37.6 


-2.7 


0.99 


5.5287 


172.1 


2 


1910+0358 


38.6 


-2.3 


2.33 


4.53 


78.8 


3 


1908+04 


39.2 


-1.4 


0.29 




217 


4 


1902+0556 


39.5 


+0.2 


0.75 


12.896 


179.7 


5 


1902+06 


39.8 


+0.3 


0.67 




530 


6 


1906+0641 


40.6 


-0.3 


0.27 


2.1352 


473 


7 


1905+0709 


40.9 


+0.1 


0.65 


4.92 


269 


8 


1901+0716 


40.6 


+1.1 


0.64 


2.40 


253 


9 


1902+07 


40.8 


+1.0 


0.49 




90 


10 


1910+07 


41.6 


-0.8 


2.71 




115 


11 


1908+07 


41.6 


-0.2 


0.21 




10 


12 


1915+07 


42.5 


-1.9 


1.54 




50 


13 


1908+0916 


43.2 


+0.4 


0.83 


0.098 


250 


14 


1904+10 


43.3 


+1.8 


1.86 




140 


15 


1854+10 


42.9 


+4.3 


0.57 




250 



measured rate with the complete system (see below), the 
five-telescope limits are quoted as the final results. 

4. Search for gamma-ray point sources 



Table 2. Known supernova remnants in the survey re- 
gion (Green 1998). Most are shell-type SNR. For W50 
individual points known as X-ray emitters are given (see 
e.g. Safi-Harb et al. 1997). For the unidentified EGRET 
source the location and la error box are given (Hartmann 
et al. 1999). 





Name 


l["] 


b[°] 


diameter [°] 


A 


3C396 


39.2 


-0.3 


~ 0.13 


B 


W50 


39.7 


-2.0 


~ 2 x 1 




SS-433 


39.71 


-2.25 






el 


39.75 


-2.66 


~0.12 




e2 


39.91 


-2.80 


~0.3 




e3 


40.09 


-3.21 


~0.45 




w2 


39.51 


-1.76 




C 




40.5 


-0.5 


~ 0.4 


D 


3C397 


41.1 


-0.3 


< 0.1 


E 




42.8 


+0.6 


0.4 


F 


W49B 


43.3 


-0.2 


< 0.1 


3EG 


3EG J1903+0550 


39.52 


-0.05 


la 0.64 
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Fig. 3. Detection rates for a point source with the flux 
and spectrum of Crab Nebula, as a function of the angle 
of incidence of the photons relative to the telescope axis. 
The simulations were performed for a five telescope sys- 
tem, and are shown for events where all telescopes trigger. 
The rates include a cut on the mean scaled width < 1.1. 
Also shown is the rate measured for the Crab Nebula, po- 
sitioned 0.5° off-axis. Full symbols: 20° zenith angle, open 
symbols: 30° 
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The cuts on the telescope images and the shower re- 
construction follow earlier work. In particular, only images 
with at least 40 photoelectrons are accepted and the cen- 
troid of the image has to be within 1.7° from the camera 
centre, in order to exclude truncated images. Showers with 
reconstructed cores up to 300 m from the center telescope 
are accepted. 

The performance of the telescope system for off-axis 
gamma-rays was investigated in detail using Monte Carlo 
simulations, for the zenith angle range in question. Fig. ^| 
shows, as a function of shower inclination relative to the 
telescope axis, the rate at which showers are triggered 
and reconstructed with all five telescopes for a Crab-like 
source. The telescope system shows a FWHM field-of-view 
of about 3.5°. Also shown is the actual Crab rate derived 
from 4.1 h of observations. During pointed observations 
with the HEGRA IACT system, the source is usually po- 
sitioned 0.5° off-center. The angular resolution is, within 
20%, independent of the inclination of the shower axis and 
varies between about 0.11° for events triggered by at least 
two of the five telescopes, to 0.07° for five-telescope events. 
The angular resolution is measured by the Gaussian width 



Fig. 4. Map of the significances. Only significances above 
1.5 sigma are shown. Superimposed are the locations of 
potential sources and the borderline of the observed region 
(see also Fig. ||) . 



of a 1-dimensional projected angular distribution, corre- 
sponding to the 40% containment radius [] Since events 
with five triggered telescopes have both a better angular 
resolution and an improved hadron rejection compared to 
events with fewer telescopes, they contribute the bulk of 
the significance in the detection of sources. In fact, indis- 
criminate combination of 2, 3, 4 and 5-telescope events 
may deteriorate the sensitivity compared to 5-telescope 
events alone. For these reasons, the current analysis is 
based exclusively on events with five triggered telescopes 

Non Gaussian tails depend on the number of telescopes 
used for reconstruction and are well below 10% for the 5- 
telescope events.. 
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Fig. 5. Upper limits for the gamma-ray flux above 1 TeV, 
in units of the Crab flux. Upper limits above the flux of 
the Crab are shown in black. 



for point sources and on four and five triggered telescopes 
for extended sources. 

For a given point source candidate, events recon- 
structed within 0.11° of the source direction were counted, 
after applying a cut on w at 1.1 to reject cosmic-ray events. 
In case of extended objects, such as SNR, the angular cut 
was increased to 0.2° or 0.25° depending on the source 
size, to ensure that the whole source is contained in the 
search region and in addition the 4-telescope events were 
used, since angular resolution is no longer critical. To es- 
timate backgrounds, three regions of the same size as the 
source region were used, rotated by 90°, 180° and 270° 
around the telescope axis, relative to the source. This 
background estimate can only be applied for sources more 
than 0.12° away from the telescope axis (otherwise the 
source region and the background regions would overlap) . 



Table 3. Upper limits for specific objects. Columns show 
the observation time used in the analysis, the number of 
ON and OFF-source events (for a three times larger OFF 
region), the significance of detection and the 99% upper 
limit in units of the Crab flux. Objects are treated as point 
sources (type P) or extended sources (type El: i? = 0.2 
deg, type E2:tf = 0.25 deg). 
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Tvne 
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ON 


OFF 




,S9% 
V 


T1 QOQ-i-OO^A 
j iyuyTU Zr04: 
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3 87 
o.o 1 


n 

u 


9 


- ± .u t 


u. iou 


T1 Q1 04-0*3^8 


p 


u.yo 


1 


a 
u 


n 7n 


n nsfi 

u.uou 


J1908-I-04 


P 


10.75 


o 


10 


-2.40 


0.046 


Ti Qn9+f)55fi 


p 


11.94 


2 


4 


0.46 


0.059 


TI 909+06 


p 


11.94 


3 


11 


-0.31 


0.060 




p 


12.46 


1 


5 


-0.49 


0.051 




p 


8.78 


2 


7 


-0.19 


0.074 


J1901+0716 


P 


7.06 


2 


5 


0.21 


0.106 


J1902+07 


P 


8.97 


2 


3 


0.74 


0.083 


J1910+07 


P 


8.40 


2 


6 


0.00 


0.087 


J1908+07 


P 


8.78 


1 


6 


-0.70 


0.068 


J1915+07 


P 


1.99 





2 


-1.07 


0.214 


J1908+0916 


P 


2.70 


2 





2.35 


0.321 


J1904+10 


P 


1.55 








0.00 


0.533 


J1854+10 


P 


1.80 





1 


-0.76 


0.285 


SNR-3C396 


P 


12.31 


4 


6 


1.04 


0.082 


SNR-G40.5 


E2 


9.62 


30 


108 


-0.90 


0.062 


SNR-3C397 


P 


8.78 


3 


4 


1.03 


0.096 


SNR-G42.8 


E2 


2.49 


12 


24 


1.12 


0.263 


SNR-W49B 


P 


2.16 


1 





1.67 


0.337 


SS433 


P 


9.02 


1 


6 


-0.70 


0.063 


SS433-el 


P 


10.00 


2 


2 


1.07 


0.080 


SS433-e2 


El 


10.00 


20 


66 


-0.38 


0.110 


SS433-e3 


E2 


7.02 


19 


60 


-0.20 


0.091 


SS433-w2 


P 


9.02 


2 


4 


0.46 


0.081 


3EG 


E2 


12.83 


39 


124 


-0.32 


0.073 



On-axis source regions are therefore excluded. A source 
region which is on-axis for one scan point will, of course, 
be off-axis for the neighbouring scan points, hence the 
net loss is small. The significance for a detection is then 
calculated according to Li & Ma (1983), with a = 1/3 
corresponding to the three background regions per source 
region. 

In a first search for gamma-ray sources at arbitrary loca- 
tions, significances were calculated for a two-dimensional 
grid with 0.0625° step size, covering the survey region. 
Note that the step size is smaller than the angular reso- 
lution to ensure that no sources are missed. This implies, 
however, that significances quoted for adjacent points are 
highly correlated. The highest significance observed is 
4.1cr, compatible with the expected distribution in the 
absence of sources. The mean significance is 0.006, the 
rms width of the distribution 1.027, which shows that the 
technique used for background estimates is sufficiently ac- 
curate. A significance map is shown in Fig. ^. The points 
with the highest significances do not correspond to lo- 
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cations of any known potential source as listed in Ta- 
bles |] and ||. Upper limits are calculated according to 
the procedure given in O. Hclene (1983) and dividing this 
number by the expected number of events for a Crab 
like source seen at the same angle 9 in the camera (see 
Fig. ||). To predict the expected event rate the HEGRA 
Monte Carlo code (Konopelko, 1999a) were used to calcu- 
late the acceptance over the FoV. In the simulations the 
Crab spectrum is taken to be d^/dE = 2.7(±0.2 ± 0.8) • 

10 -ll£-2.59(±0.06±0.1) ph cm -2 s -l Tey _i ( Konopelk0j 

1998). An additional scaling of the upper limit by 28% ac- 
counts for losses caused by the tight angular cut for point 
sources. The resulting upper limits are shown in Fig. [|. For 
the well-sampled points of the survey the limits are below 
10% of the Crab flux. The limits given assume a slope of 
the energy spectrum similar to the Crab Nebula. However, 
since limits are given for an energy which corresponds to 
the peak differential detection rate, the variation of the 
limits with the assumed spectral index is small (Aharo- 
nian et al. 1995). For a variation of ±0.5 in the spectral 
index, the limits typically change by 6%. 
Specific flux limits were derived for the potential source 
candidates, namely the 15 pulsars, the 6 supernova rem- 
nants, and the EGRET source. In case of the fairly ex- 
tended SNR W50, hot spots known from X-ray measure- 
ments were treated as point sources. The so called east- 
ern lobe (e3, see Tab. ||) and the associated knot e2 were 
treated as extended sources. W50 has been observed pre- 
viously by HEGRA and a more detailed analysis of the 
previous data set and a discussion of results on W50 will 
be given in a separate paper (Aharonian et al. 2001). 
Table |3| lists the upper limits obtained for all objects. 
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Fig. 6. Latitude dependence of the diffuse emission mea- 
sured by EGRET for photon energies above 1 GeV, in the 
range of Galactic longitude 30° to 50°, together with a 
fit by a sum of two Gaussians with widths and centres of 
0.95°, 3.16° and —0.05°, 1.5°, respectively, and an inten- 
sity ratio of 4:1 in amplitude. A value of 0.1 was taken for 
the constant background. Data points taken from Hunter 
et al. 1997. In addition, Gaussian profiles with a a of 0.5° 
,1.0°, 2.0° and 3.0° are shown. 



5. Search for diffuse gamma-ray emission from 
the Galactic plane 

Compared to the search for point sources, the search for 
diffuse gamma-ray emission from the Galactic plane is 
complicated by the extended structure of the emission re- 
gion. The structure may be extended in latitude beyond 
the field of view. Fig. ^ illustrates the profile in Galactic 
latitude as measured by EGRET in the longitude range 
30° to 50°, at energies above 1 GeV. The EGRET lat- 
itude profile can be represented as a sum of two Gaus- 
sians and an constant background value. Since the angu- 
lar resolution of the EGRET instrument above 1 GeV is 
narrower than the latitude extend of the diffuse emission 
the latitude dependence in Fig. ^| indicates real structure. 
Because of the possible change of the primary emission 
mechanism between the GeV and TeV range, from 7r° de- 
cay to inverse Compton scattering (Porter & Protheroe, 
1997; Pohl & Esposito, 1998), and because of the com- 
pletely different parent and target populations in the two 
cases, the latitude dependence may differ at TeV energies. 
At GeV energies, the latitude distribution of gamma-ray 
emissivity is characterized by the scale height of inter- 



stellar gas, around 100 pc for atomic hydrogen (Lockman 
1984) and somewhat less for molecular hydrogen (Dame 
1987). In contrast, inverse-Compton interactions with the 
uniform microwave background photons are governed by 
the scale height of the electron component of cosmic rays, 
which might be characterized by the kpc scale describing 
the coupling of cosmic rays to matter. The scale height of 
far-infrared target photons, on the other hand, was given 
by Cox, Kriigel and Metzger (1986) as 100 pc, similar to 
the scale height of the gas. Porter and Protheroe (1997) 
find, from numerical simulations of electron propagation, 
a scale height of more than 500 pc for the inverse-Compton 
emissivity at 1 TeV. 

In the following, we will discuss three techniques to search 
for diffuse gamma-ray emission. The techniques differ in 
the degree of assumptions they make concerning the lati- 
tude dependency of the diffuse emission, and also in their 
sensitivity to systematic variations in the performance and 
characteristics of the telescopes. 

The most robust and model-independent - but also least 
sensitive - technique to derive limits on the diffuse flux 
simply selects events according to their shapes as gamma- 
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Fig. 7. Difference between the distributions in mean 
scaled width for the on-region (|6| < 2°) and for the off- 
region > 2°). The dashed lines indicate the expected 
gamma-ray region. 



ray candidates. A cut on w less than 1.0 keeps about 
1/2 of the gamma-rays, but rejects cosmic rays very effi- 
ciently (see Fig. ^). In order to achieve the best separation 
between gamma-ray images and cosmic-ray images, only 
five-telescope events were used. Such events with small w 
include genuine gamma-rays, electron showers from the 
diffuse cosmic-ray electron flux, and the tail of the dis- 
tribution of cosmic-ray nuclear showers. Assigning all 428 
events after cut as diffuse gamma-rays a 99% upper limit 
on the diffuse gamma-ray flux at 1 TeV of 23.4 -10 -15 ph 
cm~ 2 s _1 sr _1 MeV -1 results, for |6| < 5°. As mentioned 
in the discussion of the point source limits, the measure- 
ment determines the integral flux above the energy thresh- 
old of the telescope system, rather than directly determin- 
ing the differential flux. Therefore the result depends in 
principle on the assumed spectral index. However, since 
flux values are quoted at energies corresponding roughly 
to the peak detection rates, the limits vary only very little 
with the spectral index. We note that the diffuse electron 
flux (see, e.g., the compilation by Barwick et al. 1998) 
should contribute about 1/3 of the number of gamma- 
ray candidates; the limit can of course also be viewed 
as a limit on the electron flux, since the selection cuts 
are equally efficient for gamma-ray induced showers and 
electron-induced showers. 

The limit obtained by this technique can be improved by 
subtracting, on a statistical basis, the number of gamma- 



ray candidate events resulting from cosmic-ray electrons 
or protons. Since it is virtually impossible to verify that 
simulations properly account for the tail towards very 
small w(< 1.1) of the distribution of proton showers, such 
a subtraction has to be performed using an experimental 
background region, sufficiently far away from the Galac- 
tic plane such that diffuse gamma-ray emission from the 
plane is most probably negligible. Such a background data 
sample will only contain the isotropic electron and nu- 
cleoli flux. In order to minimize instrumental effects, the 
background sample should be taken at the same time, 
and at identical zenith angles. Data sets can be normal- 
ized to each other using the rates of events with large 
w (> 1.4), well outside the gamma-ray region. Unfortu- 
nately, the availability of suitable background data sam- 
ples with the same telescope configuration is limited to 4.1 
h of data, and the statistical error of the background data 
set limits the sensitivity. After a w < 1.1 cut 1928 events 
survived compared to 141 events in the referenc sample. 
The scaling factor between the two datasets is determined 
to 13.3 and the MC simulations predict 555 events for a 
Crab like source smeared out over a FoV of 1.5°. After 
subtraction of isotropic components, we find a 99% limit 
on the diffuse flux in the region \b\ < 5° of 10.4 -10 -15 ph 
cur 2 s" 1 sr" 1 MeV" 1 at 1 TeV. 

The final, and most sensitive analysis makes the assump- 
tion that diffuse gamma-ray emission from the Galactic 
plane is limited to the central parts of the scan region, 
and uses the outer parts of the scan region to estimate 
backgrounds. The range \b\ < 2° was considered the sig- 
nal region, the range |6| > 2° the background region. This 
cut is close to optimal for emission profiles with an rms 
width between 1° and 2° and results in a good balance 
in observation time (awl) between the ON and OFF 
regions. To ensure optimum quality of the events, only 
four- and five-telescope events were used, and the field of 
view was restricted to 1.5° from the optical axis. A cut 
at 1.1 on the w was applied to reject cosmic-ray back- 
ground 0. To account for a possible zenith-angle depen- 
dence of background rates, data were grouped into four 
different ranges in zenith angle, 20°-24°, 24°-28°, 28°-32° 
and 32°-36°. Also, the analysis was carried out separately 
for each scan band. For each range in zenith angle and 
each scan band, the expected number of events in the sig- 
nal region was calculated based on the number of events 
observed in the corresponding areas of the camera for the 
background region. The expected and observed numbers 
of events were then added up for all zenith angles and scan 
bands. With a total number of 2387 gamma-ray events in 
the signal region, compared to 2353 expected events, there 



2 The cuts on w differ between this and the first analysis; 
here, the goal is to optimize the significance of a weak signal 
S, e.g., the ratio S/yB, whereas analysis of the first type, 
where all backgrounds B are counted as potential signal, need 
to optimize S/B. 
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is no significant excess. As a cross check, Fig. [T] shows the 
background-subtracted distribution in w for events in the 
|6| < 2° region. There is no indication of a significant ex- 
cess in the gamma-ray region around a w of 1 indicated 
by the dashed lines. Also for larger values of w, there is 
no significant excess or deficit, showing that background 
subtraction works properly. 

In order to translate the limit in the number of events 
into a flux limit, one now has to make assumptions con- 
cerning the distribution in Galactic latitude of the diffuse 
radiation, since a spill-over of diffuse gamma-rays into the 
background region |6| > 2° will effectively reduce the sig- 
nal. For a profile with a width less or comparable to the 
EGRET profile - about 1° rms - a correction of 12% is 
applied and one finds a limit 6.1 • 10~ 15 ph cm -2 s -1 sr -1 
MeV -1 for the diffuse gamma-ray flux at 1 TeV, averaged 
over the \b\ < 2° region. The limit refers to an assumed 
spectral index of -2.6, and changes by +13% for an index 
of -2, and by -5% for an index of -3. For wider distribu- 
tions of 2° and 3° rms, the limit changes to 8.2 ■ 10 -15 ph 
cm~ 2 s" 1 sr- 1 MeV" 1 and 12.8 • 1CT 15 ph cm~ 2 s -1 sr -1 
MeV -1 , respectively. Fig. |^ compares the upper limit with 
the extrapolation of the EGRET flux; also included are the 
limits by the Whipple group (Reynolds et al. 1993; LeBo- 
hec et al. 2000) and of the Tibet array (Amenomori, 1997). 
Connecting the EGRET points with the HEGRA upper 
limit (and ignoring the highest-energy EGRET point with 
its large errors) one finds a lower limit of 2.5 on the spec- 
tral index of the diffuse gamma-ray emission. 

Model predictions attempting to explain the excess 
flux in the GeV region by assuming an increased inverse 
Compton component (Porter & Protheroe 1997, Pohl & 
Esposito 1998), or by taking the contribution from unre- 
solved SNRs into account (Berezhko & Volk, 2000), are 
given in the literature for different ranges in Galactic lati- 
tude and cannot be directly compared without assum- 
ing a latitude dependence of the diffuse flux. In Fig. || the 
solid line gives the model prediction of Berezhko & Volk 
2000 scaled by a factor of 3, taking into account the differ- 
ent latitude and longitude ranges used for the upper limit 
and the model and using the EGRET measurements at 20 
GeV as a guideline. 

Even for pessimistic assumptions of a rather wide lat- 
itude dependence the model of Pohl & Esposito (1998), if 
the inverse Compton flux is extrapolated from the 50 GeV 
range discussed in the paper to the TeV range, exceeds the 
HEGRA limit. At TeV energies, one is still sufficiently far 
from the Klein-Nishina regime such that the power-law 
extrapolation should be valid (Porter & Protheroe 1997). 
Of course, a break in the electron injection spectrum could 
cause a corresponding break in the gamma spectrum be- 
tween 50 GeV and 1 TeV, and could be used to make the 
model consistent with the experimental limit. 
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Fig. 8. Upper limits for the diffuse gamma-ray flux de- 
rived by this experiment. Assuming that all detected 
events are gamma-rays (1), using a independent data set 
for background subtraction (2), and using \b\ > 2° as back- 
ground data and assuming the spatial distribution mea- 
sured by the ERGRET instrument (38° < I < 43°, |6| < 
2°) (3). For detailed explanations see text. Also shown is 
the EGRET flux for 35° < I < 45°, \b\ < 2°, the Whip- 
ple upper limits for 38.5 < I < 41.5°, |6| < 2° (Reynolds 
et al. 1993, LeBohec 2000) and the Tibet upper limit 
(Amenomori, 1997). The dotted line indicates an extrap- 
olation of the EGRET data with an index of 2.5. The 
dashed line indicates the scaled 'leaky box' model predic- 
tion by Berezhko & Volk, 2000 (BV2000) for the spatially 
averaged gamma-ray emission from the Galactic plane by 
those Galactic CRs that are still confined in their unre- 
solved SNR sources, scaled to |6| < 2°. 



6. Concluding remarks 

A survey of the region -10° < b < 5°, 37° < / < 43° near 
the Galactic plane did not yield evidence for TeV gamma- 
ray point sources, with typical flux limits of 10% of the 
Crab flux. In particular, 15 pulsars, 6 supernova remnants 
and one unidentified EGRET source were not detected as 
strong TeV sources. 

A search for diffuse gamma-ray emission resulted in 
an upper limit of 6.1 • 10~ 15 ph cm -2 s -1 sr -1 MeV -1 
at 1 TeV, averaged over the region 38° < I < 43°, |6| < 
2° and assuming the spatial emission profile measured by 
the EGRET instrument. Since the analysis used to derive 
this limit is only sensitive to the variation of the diffuse 
flux with 6, rather than its absolute value, a distribution 
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significantly wider than at EGRET energies will increase 
the limit. Other variants of the data analysis are sensitive 
to the absolute flux, but give less stringent limits of 23.4 
and 10.4-1CT 15 ph cm" 2 s" 1 sr" 1 MeV" 1 . The limit on the 
TeV gamma-ray flux can be used to derive a lower limit 
on the spectral index of the diffuse radiation of 2.5, and 
to exclude models which predict a strong enhancement 
of the diffuse flux compared to conventional mechanisms. 
However the TeV flux limit is only about a factor of f .5 
larger than the predicted flux from unresolved SNRs. A 
more sensitive survey should therefore be able to test this 
prediction, together with a determination of the "diffuse" 
TeV gamma-ray spectrum that is directly related to the 
Galactic CR source spectrum. 
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